Spin excitations in the excitonic spin-density-wave state of the iron pnictides 
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Motivated by the iron pnictides, we examine the spin excitations in an itinerant antiferromagnet 
where a spin-density wave (SDW) originates from an excitonic instabihty of nested electron-hke 
and hole- like Fermi pockets. Using the random phase approximation (RPA), we derive the Dyson 
equation for the transverse susceptibility in the excitonic SDW state. The Dyson equation is solved 
for two different two-band models, describing an antiferromagnetic insulator and metal, respectively. 
We determine the collective spin-wave dispersions and also consider the single-particle continua. The 
results for the excitonic models are compared with each other and also contrasted with the well- 
known SDW state of the Ifubbard model. Despite the qualitatively different SDW states in the two 
excitonic models, their magnetic response shows many similarities. We conclude with a discussion 
of the relevance of the excitonic SDW scenario to the iron pnictides. 

PACS numbers: 75.30.Fv, 75.10.Lp 



I. INTRODUCTION 

The recent discovery of superconductivity in iron pnic- 
tides has sparked a tremendous research effort iiS. The 
remarkably high superconducting transition temperature 
Tc of some of these compounds,'^ their layered quasi-two- 
dimensional structure,^ the proximity of superconductiv- 
ity and antiferromagnetism in their phase diagramsj^i^'^ 
and the likely unconventional superconducting pairing 
state^ i^^d^d^d'^ are reminiscent of the cuprates.^"* It is 
a tantalizing prospect that the iron pnictides can shed 
new light onto the problem of unconventional high-Tc 
superconductivity in general. 

For this it is essential to assess the differences be- 
tween the cuprates and the iron pnictides. For exam- 
ple, the pnictides have a much more complicated Fermi 
surfaceJ^ The antiferromagnetic states in the two fami- 
lies are also qualitatively different. In the cuprates, su- 
perconductivity appears by doping an insulating antifer- 
romagnetic parent compound. The pnictide parent com- 
pounds iJFeAsO {R is a rare-earth ion) and j4Fe2As2 {A 
is an alkaline-earth ion) are also antiferromagnets, but 
there is compelling evidence that they display a metallic 
SDW state: the value of the magnetic moment at the Fe 
sites is smallr^ii^ the compounds display metallic trans- 
port properties below the Neel temperature T/vr^'^'^ 
and ARPES and quantum oscillation experiments find a 
reconstructed Fermi surface below TatJ^iSS 

The electron-phonon interaction in the pnictides is 
much too weak to account for the high Tc values. 
Instead, the most likely candidate for the "glue" 
binding the electrons into Cooper pairs are spin 
fluctuationSfiiii^iii^ which are enhanced by the proxim- 
ity to the SDW state. A proper understanding of the 
SDW phase is therefore likely the key to the physics 
of the pnictides. Intriguingly, ab-initio calculations sug- 
gest that the nesting of electron and hole Fermi pockets 
is responsible for the SDW,^^ indicating that, like the 
superconductivity, the antiferromagnetism of these com- 
pounds has a multiband character. The best known ma- 



terial where a SDW arises from such a nesting property is 
chromiu m^^'^'^1^^1^^1^^ and this mechanism has also been 
implicated for manganese alloys 3i 

The SDW in these compounds belongs to a broader 
class of density-wave states. Consider a material with 
electron-like and hole-like Fermi pockets separated by a 
nesting vector Q in the presence of interband Coulomb 
repulsion. Performing a particle-hole transformation on 
one of the bands, we obtain an attractive interaction be- 
tween the particles in one band and the holes in the 
other. Within a BCS-type mean-field theory, the at- 
tractive interaction causes the condensation of interband 
electron-hole pairs (excitons) with relative wave vector 
Q, thereby opening a gap in the single-particle excita- 
tion spectrum. ■^^ Although the interband Coulomb repul- 
sion causes the excitonic instability, additional interband 
scattering terms are required to stabilize one of several 
different density- wave states, such as a SDW or a charge- 
density wave (CDW),^i20|3i 

Several authors have discussed the SDW state of the 
pnictides in terms of an excitonic instability of nested 
electron and hole Fermi pockets without regard to the or- 
bital origin of these hands M^^^i^^j^^i^^^ An alternative 
school of thought emphasizes the importance of the com- 
plicated mixing of the iron 3c? orbitals at the Fermi energy 
and of the various inter-orbital interactionsi^'^i^'^ 
These two approaches are not contradictory, however, 
since the excitonic model can be understood as an ef- 
fective low-energy theory for the orbital modelsi^i^ Fur- 
thermore, even in an orbital model, the SDW state is still 
driven by the nesting of electron and hole Fermi pockets. 
Indeed, at the mean-field level all these models yield qual- 
itatively identical conclusions. A conceptually different 
picture based on the ordering of localized moments has 
also been proposedJ^i^i^i^ Although it is hard to rec- 
oncile with the observed metallic propertiesi^iii and the 
moderate interaction strengthsf^^i^ this picture is consis- 
tent with several neutron-scattering experiments.— At 
present, it is difficult to discriminate between the itiner- 
ant (excitonic) and localized scenarios, as the dynamical 
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spin response of the itinerant models is unknown. It is 
therefore desirable to determine the spin excitations in 
the excitonic SDW model. 

It is the purpose of this paper to examine the trans- 
verse spin susceptibility within the excitonic SDW state 
of a general two-band model. We work within the lim- 
its of weak to moderate correlation strength, using the 
RPA to construct the Dyson equation for the suscepti- 
bility. In order to understand the generic features of the 
spin excitations in the excitonic SDW state, we calculate 
the RPA susceptibilities for the simplest model showing 
this instability. We pay particular attention to the spin 
waves (magnons) and damped paramagnons. In the sim- 
plest model, however, the SDW state is insulating. We 
therefore verify the robustness of our results by applying 
our theory to a system where portions of the Fermi sur- 
face remain ungapped in the SDW phase, as in the iron 
pnictides. We contrast our results for the excitonic SDW 
state with those for the SDW phase of the single-band 
Hubbard model, which is commonly used to describe the 
antiferromagnetic state of the cuprates. 

The structure of this paper is as follows. We com- 
mence in Sec. |lT] with a brief review of the RPA-level 
results for the transverse susceptibility in the SDW state 
of the Hubbard model4^i^i^ We then proceed in Sec. IIIII 
with a general discussion of the excitonic SDW state in a 
two-band model and present the Dyson equation for the 
transverse susceptibilities. The RPA susceptibility and 
spin-wave dispersion is then calculated for the insulating 
and metallic excitonic SDW models in Sees. IIII Al and 
IIIIBl respectively. All presented results are calculated in 
the limit of zero temperature. In order to properly com- 
pare the different models, we choose interaction strengths 
such that the zero-temperature SDW gap is the same. We 
conclude with a comparison with experimental results in 
Sec. II VI and a summary of our work in Sec. |Vl 



II. HUBBARD MODEL 



The Hamiltonian of the Hubbard model reads 



k,o- 



~ 5Z ^kCk g-Ck^g. -I- — ^ C^q^t'^k.tCk'-qa'^k'a: (1) 
k,k',q 



where cj^. ^ {c^ ^) creates (destroys) an electron with mo- 
mentum k and spin a. We assume a two-dimensional 
(2D) nearest-neighbor tight-binding dispersion = 
—2t {cos k^a -\- cosfcyo) where a is the lattice constant. 
We plot the band structure and the resulting Fermi 
surface at half filling in Fig. [TJa) and (b), respectively. 

At half filling and sufficiently low temperature T, the 
Hubbard model is unstable towards a SDW state with 
nesting vector Q = {n/a^n/a), which connects opposite 
sides of the Fermi surface. We assume a SDW polar- 
ized along the z axis, and decouple the interaction term 
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FIG. 1: (Color online) (a) Band structure and (b) Fermi sur- 
face of the Hubbard model for (7 = 0. In (b), the nesting 
vector Q = {ji/a,-K/a) is also shown. 



in Eq. (jlj by introducing the SDW gap 
Y^' / t 

' / . ^ \ '^k+Q,o-Ck,o- 



V' 



(2) 



k,(T 



The primed sum denotes summation only over the re- 
duced, magnetic Brillouin zone. Diagonalizing the mean- 
field Hamiltonian, we find two bands in the reduced Bril- 
louin zone with energies -E±,k = ±-\/e| + A^. In the 
following, we will assume t = 1 eV and U = 0.738 eV, 
which gives a critical temperature for the SDW state of 
TsDw = 138 K and a T = gap A = 21.3 meV. 
The dynamical spin susceptibility is defined by 

X,,(q,q';*c.„) = ^ J^dT {TrS\q,r)S^{-q',0))e'^"^, 

(3) 

where is the time-ordering operator and 

- E E % ^k,.' (r) ■ (4) 

k s.s' 

Because of the doubling of the unit cell in the SDW 
state, the susceptibility in Eq. ^ is non-zero for q = q' 
and q = q' -f Q, the latter referred to as the umklapp 
susceptibility.'*® Both appear in the ladder diagrams for 
the transverse susceptibility, yielding the Dyson equation 

x-+(q, q';«^n) 
= 6q,q'X^-l{^, q; it^n) + ^q+Q,q'X-°l (q, q + Q; it^n) 

+ C/xL°|(q, q; iw„)x_+(q, q'; it^n) 

;7xL°|(q,q+ Q;it^„)x_+(q + Q,q';iw„), (5) 

where the superscript (0) indicates the mean-field sus- 

ceptibilties. Explicit expressions for y l(c|, q; iLUn) and 

X h(q) q + Q; *'^n) can be found in Ref. |50- 

We plot the imaginary part of x i-(qi^) = 

X-+(q>q;w) along the line q = iqx,qy = Qx) in Fig. H 
The calculation of the mean-field susceptibilities in the 
Dyson equation (O was performed over a 10000x10000 
k-point mesh. In the analytic continuation icun to + i6 
we assume a finite width (5 = 1 meV. Smaller values of 6 
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FIG. 2: (Color online) Imaginary part of the transverse susceptibility in the Hubbard model for q = {<lx,<ly ~ Qx)- The 
spin- wave dispersion is visible as the dark line running across the figure at a; < 2 A = 42 meV. Note the logarithmic color scale. 



and finer k-point meshes do not produce qualitative or 
significant quantitative changes in our results. 

Imx in Fig. [5] displays very different behavior 

for energies oj < 2 A = 42.6 meV and ut > 2 A. In the for- 
mer region, the dispersion of the collective spin waves is 
clearly visible as the sharp dark line. The finite width of 
this hne is a consequence of the broadening S. The disper- 
sion is almost flat for 0.1 n/a ^ Qx — Qy ^ 0.9 7r/a, where 
it lies very close to uj = 2A. The distribution of spec- 
tral weight for the spin wave is asymmetric, with much 
greater weight close to q = Q than at q = 0, reflecting 
the suppression of long- wavelength spin excitations in the 
SDW state,^ 

For UJ > 2A, we find a continuum of excitations. It 
starts abruptly at w = 2A, corresponding to the min- 
imum energy for a single-particle excitation across the 
SDW gap. This minimum is the same at all k points 
lying on the Fermi surface shown in Fig. [Ub). By in- 
spection, we see that for every value of q, there exist 
points k and k -I- q lying on the Fermi surface so that the 
minimum energy required for any excitation is lu = 2A. 
We also see that Im x ^) tends to decrease with in- 
creasing Lj. This can be understood in terms of the den- 
sity of states (DOS) in the non-interacting model: the 
DOS has a van Hove singularity at the Fermi energy and 
decreases monotonically as one moves to higher or lower 
energies. For an occupied state with energy ujo below the 
Fermi energy, the density of unoccupied states with en- 
ergy uju above the Fermi energy therefore decreases with 
increasing lu — Uu — lUq, and hence the "density of exci- 
tations" contributing to the transverse susceptibility also 
decreases with increasing lu. 

Close to q = 0, the continuum is bounded from above 
by the line lu — -q where is the Fermi velocity along 

k = {kx,ky = kx). The peak in Imx i-(Qi^) at this 

edge of the continuum is due to single-particle excitations 
across the Fermi energy in the same branch of the band 
structure. A rather weak dispersing feature also appears 
within the continuum near q = Q, as shown in more 
detail in Fig. [3l This paramagnon originates from single- 
particle excitations into the back-folded band. Like the 
feature at small q, the paramagnon disperses with the 
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FIG. 3: (Color online) Imaginary part of the transverse spin 
susceptibility in the Hubbard model for q = {qx,qy = qx) 
close to Q. Note the linear color scale. 



Fermi velocity. The paramagnon and spin-wave disper- 
sions curve away from one another in what appears to be 
an avoided crossing. 

Solving Eq. (O for x ^{'ii^^'i^n) requires the inver- 
sion of a 2 X 2 matrix. The determinant I?(q, ioun ) of this 
matrix is 



1 - c^x-°l(q, q; i^n) 

l-[/xL"i(q + Q,q+Q;ic^„) 



(6) 



Making the analytic continuation iLUn — > w -t- iO+, the 
solution of Re2?(q, w) = yields the spin- wave disper- 
sion. At low energies, it has a linear dependence upon 
(5q = Q— q, i.e., uj = csw li^ql, where csw is the spin- wave 
velocity. An expression for cgw is obtained by expand- 
ing I?(q, io) about q = Q and lu = ,48,49,50 jj^ agreement 
with Ref. Ea, we find 



csw 



'-4(l/[/- A2x)t27 



x/U 



(7) 
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FIG. 4: (a) Spin-wave velocity csw in the Hubbard model 
as a function of U for t = 1 eV. (b) Comparison of the spin- 
wave dispersion and the low-energy linear form as a function 
of 5q = Q - q for 17 = 0.738 eV. 



where 



1 1 
^? in 



(8) 



(cos k^a + cos kxtt cos kyo) 



—7 sm krtt 



(9) 



The spin-wave velocity is plotted as a function of U 
in Fig.m^a), while we compare the low-energy linearized 
form of the spin-wave dispersion to the numerically- 
determined result in Fig. IH^b). As can be seen, the lin- 
earized result holds only for small energies lo < 0.5 A. 



III. EXCITONIC MODEL 

In this section we discuss the excitonic SDW in a gen- 
eral two-band model with Fermi-surface nesting. We 
begin by outlining the known results for the mean-field 
SDW state..2i>25i28^^ ^^.j^e Hamiltonian as 



H — Hq + Hi, 
where the non-interacting system is described by 



(10) 



Hn = 



E E [(^k - ^^)cLc^. + (4 - M)/L/k.] (11) 



interband Coulomb repulsion, 



H, 



cf 



9ci 
V 



E E "^k+q.o-Ck.CT/k'-q.cr'/k'.CT'' (1^) 
k,k' ,q a",cj' 

and two distinct types of correlated interband transitions, 



^iTa- ^ E (4+q,T4'-q,i/k'.i/k,T+H.c.), (15) 

k,k',q 

HlTh= ^ E E"=k+q,^/k'-q,^'Ck%<T'/k,<T- (16) 
k,k' ,q (7,(t' 

The interband interaction terms are responsible for a 
density- wave instability when the electron and hole Fermi 
pockets are sufficiently close to nesting. A number of dif- 
ferent density- wave states are possible*^ a CDW with ef- 
fective coupling constant gcBW — Sc/— 52a — 2 (72b, a SDW 
with coupling gsDW = 5c/ + 52a, a charge-current-density 
wave (CCDW) with 5ccdw — 5c/+52a— 252b, and a spin- 
current-density wave (SCDW) with 5scdw = 5c/ — 52a- 
In order to model the iron pnictides, we henceforth as- 
sume that the SDW state has the largest coupling con- 
stant. 

In the presence of a SDW polarized along the z- 
axis and with nesting vector Q, the effective mean-field 
Hamiltonian is written as 

i/MF = i?0 + E E (4,./k+Q,. + H.C.) , (17) 
k cr 

where the excitonic gap 

» 5SDW 



2V 



EE^(^k,-^k+Q 



(18) 



k cr 



is assumed to be real. The precise relationship between 
A and the magnetization is somewhat complicated i^^i'^° 
To elucidate it, we define the field operator. 



(r) = E ['^kx(r)ck,^ + ¥'k,/(r)/k,^] 



(19) 



where ip]i^a(j) is a Bloch function for the band a. The 
local magnetization M(r) is then 



and cjj. ^ (/k cr) creates an electron with spin a and 
momentum k in the electron-like c band (hole-like / 
band). The second term in Eq. pH)) describes the in- 
teractions in the model system. Following Ref.s [3l| 
and [33, we take this to consist of five on-site terms 
Hi — Hcc+H ff+Hcf+HiTii+HiTh that arise naturally in 
the low-energy effective theory of a multi-orbital model. 
These correspond to intraband Coulomb repulsion. 



M(r) = -^EE E ^k,a(r)^k^.(r) 



s.s' k.k' a,b—c,f 



X e 



-i(k-k': 



t CTs.s' , 



(20) 



where 5 is the (7-factor and [Ib is the Bohr magneton. 
Only in the limit when </5k,Q(r) is constant do we find A 
to be simply related to the magnetization. 



Hnn — 



V 



E 

k,k',q 



'^k+q,T'^k,Tcj(.,_q |Ci(./_|, 



(12) 



Hff - /.l+q.T/k,T/k'-q,i/k'.i' (13) 



k,k',q 



M(r) — (cos Q • r) Bz 

5SDW 



(21) 



For simplicity, we follow Refs. Il3ll23ll27ll52l in assuming 
constant Bloch functions. 
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In calculating the susceptibilities, we make use of the 
single-particle Green's functions of the mean-field SDW 
state. The two normal (diagonal in band indices) Green's 
functions are defined by 



,(22) 



k+Q 



(23) 



while the anomalous (band- mixing) Green's functions are 

r0 



aA 



<.(*^n) = rf^(r.Ck,.(r)/^+Q^,(0))e^' 



= G{+q (iw„). 



(24) 



(25) 



The functions E±,]^ are the dispersion relations for the 
reconstructed bands, 



E. 



1 



±,k 



fi' 



(26) 



For energies much larger than A we have -E+,k ~ ^k+q 

and -E-,k ~ e^. 

The total spin operator is written as 

S(r) = ^v4(r)^V^..(r) 



a,b—c,f k.q s.s' 
^ ^ a,b=cj q 



(27) 



where Sa,6(q) is a generalized spin operator. The dy- 
namical spin susceptibility is then defined by 

Xji(q, q';iuj.n) 
= fEE r(r.5:,,(q,r)5^,,,(-q',0) 

= Y.Y.xt''''(.^,<l';^^n). (28) 
a, 6 a' ,b' 

The generalized susceptibilities xfj"^ ^ (q, q'; *i^n) are cal- 
culated using the RPA. We are only concerned with the 



transverse susceptibility, which is obtained by summing 
the ladder diagrams. This yields the Dyson equation 



aba'b' 
A — h,00 



i- / c c abba (0) , r- r- abba (0) 

= dci,ci' [Oa',bOb',aX^ + ^ 00 + Oa',A.aX- + , 00 

, e / 1: r abba (0) , c- c abba (0)\ 

+ dq+Q.q' [da',bdb',aX^ + ^QCl + a' .b'^b' ,aX- + , QQ ) 

_1_ ^ ( ^,"-bcc{0) cca'b' , abcc{0) cca'b' ^ 
-r i/cc yX — ^,00 A — 1-,00 ^ A — h.OQA — h,QO^ 



„ {abff{0) ffa'b' ,abff{Q)ffa'b' 

' yff yX — h,oo A — h^oo'T-A — h,oQA — h,Qo 

f abcfiO) fca'b' ,abcf{0) fca'b' 
' 9cf \ X- + , 00 A- + , 00 + A- + , OQ A- + , QO 

abfciO) cfa'b' abfc{a) cfa'b' \ 

A- + ,00 A- + ,00 ^ A- + ,0QA_ + ,Q0 J 

f abcfiO) cfa'b' , abcfiO) cfa'b' 
- '/2a yX — h, 00 A — h, 00 "I" A — OQ A — ^, QO 

q6/c(0) /ca'6' abfc(O) fca'b' \ 

A — ^,00 A — h,00 ^ A — h,OQA — h,QO J 

f abcc{0) ffa'b' abcc{0) ffa'b' 

■ 92h \ X^ + , 00 X- + , 00 + X_ + ^ 0QX- + , QO 

, abff (a) cca'b' . abff (a) cca'b' \ (cyQ\ 

+ X-+,oo X-+,oo + X-+,oqX-+,qo ) ' l^yj 



where we have adopted the short-hand notation 

aba'b' (0) aba'b' (Q)/ ■ <. /on\ 

X-+,mn ^ X-+ ' '(q + m,q + n;itj„), (30) 



X 



aba'b' 
- + , inn 



x!,Y'''(q + m,q' + n;jc^„). (31) 



Note that X-+,mn'' does not depend on q'. We have 
also introduced the notation a = c{f) when a = f {c). 
The first two lines of Eq. give the mean-field suscep- 
tibilties, obtained by using Wick's theorem to contract 
the correlation function in Eq. ((28|) into products of two 
mean-field Green's functions. On the first line of Eq. 
we find the correlators resulting from the product of two 
normal Green's functions, 



abba (0) 

X — h,oo 



Y E ^ E GL';T(*^n)Gr+q,i(«^« - *^»). 

k ' w„ 



and from the product of two anomalous Green's func- 
tions. 



abba (0) 



X-+,oo' = E ^ E G^'T(^^")Gr+q,i(*^» - ^^n). 

k ^ w„ 



On the second line we find the umklapp susceptibilties 
which, as in the Hubbard model, are the product of a 
normal and an anomalous Green's function: 



x-'S = )^^)^G^':T(*^«)Gr+qa(*^n-^^n). 

k iVn 

x-SoS = -7 E ^ E G^'T(*-")Gr+qa(-n - 

k w„ 

The remaining lines of Eq. (j29p give the ladder sums for 
the various interactions: on the third and fourth lines 



1 \ ' 1 \ ' r'bb 
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we have the intraband Coulomb interactions, on the fifth 
and sixth lines the interband Coulomb interaction, and 
on the last four lines the two types of correlated transi- 
tions. In Fig.[n{a) and (b), we show a diagrammatic rep- 
resentation of the Dyson equation for x'^-i-'^oo ^^'^ X'^+oo: 
respectively. Note that the Dyson equation is also valid 



abba (0) 



are the 



in the normal state, in which case the x ^ oo 

only non-zero mean-field susceptibilities. 

Prom the structure of the Dyson equation, we see that 
X°^+ mn is only non-zero for q = q' and m, n G {0, Q}. 



We observe that Eq. (|29|) may then be written as four 
independent sets of coupled equations for 

/v'^'^'^^ v^^'"' -v"^'"'' \ C^^a) 

— h,OOJ A — h,00' A — |-,Q0' A — h,QOj' yOAA) 

f //// ccff cfff fcff \ ronu\ 
|X-+.oo> X-+,oo> Qo- X-+,Qo/ ' (32b) 

{fccf Cfcf CCCf ffcf \ for) \ 

X — 1-,00' X — h,00' X — h,QO' X — h,QOJ ' [OZC) 

{^ff^ f^f^ ccfc fffc 1 f'QOrl^ 

X-+,00' X — 1-,00' X^+^QQi X_-|-^Qoj ■ (.iza) 

Note that this includes x"^+ mn fo'' — QQ ^i^'i 

mn = OQ by symmetry; all other possible transverse 
susceptibilities vanish. The first two sets contain the con- 
tributions to the intraband susceptibility, which involve 
spin-flip transitions within the c and / bands: 



^ bf 

a.b—c.j 



J^^ dT (TrS-M, t)5+ (-q, 0)) e-"- 



— X-+,oo + X-+,oo + X-+,oo + X-+,oo- y-^"^) 

The last two sets contain the contributions to the in- 
terband susceptibility, which involve spin-flip transitions 
between the c and / bands: 



x'LT(q,*^n) 

V ^ 

a,b—c,f 



f dr(TrS--Sci,r)S+-^i-ci,0) 

J 



- ^C//C cfcf fcfc fccf ,oA\ 

We note that the Dyson equation for the interband sus- 
ceptibilities has been previously obtained in Refs. [13 
and [52 for the case where only the interband Coulomb 
interaction Eq. (fT4|) is non-zero. From Eq. (|28|) we see 
that the total transverse susceptibility is the sum of the 
intraband and interband contributions, 

x_+(q, liOn) = x'lT(q> ^^n) + x'LTXq, *^«)- (35) 

Since the interband and intraband susceptibilities involve 
qualitatively different types of excitations, considering 
these seperately offers greater physical insight into the 
magnetic response than the total susceptibility. 

In the following sections we discuss the transverse sus- 
ceptibility for two different models of the band struc- 
ture. For simplicity, we restrict ourselves to the case 



9cc = 9ff = 92h = 0, as these interactions do not drive 
the SDW instability. We emphasize, however, that the 
preceeding results are valid for any choice of couplings in 
both the normal and SDW states. Except where stated 
otherwise, we furthermore set g2a = 0, as at reasonable 
coupling strengths we find very little change in the trans- 
verse susceptibility upon varying gcf and g2a while keep- 
ing gsDW = gcf + 52a fixed. Unless explicitly mentioned, 
we have used a 10000x10000 k-point mesh and a width 
(5 = 1 meV to calculate the mean-field susceptibilities. 



A. Insulating SDW state 



We first examine an excitonic model with perfect nest- 
ing between the electron and hole bands, i.e., ej^ = — Ek+q 
for all k. Although hardly realistic, at the mean-field 
level it exactly maps onto the BCS model after particle- 
hole transformation.^- It is therefore useful for obtain- 
ing physically transparent results and is frequently en- 
countered in the literature i^^'^Sii^'^i^ We assume the 
2D band structure 



ej^ — 2t {cos k^a + cos ky a) + Co, 
— 2t {cos k^a + cos kyu) — eo, 



(36a) 
(36b) 



where we set t = 1 eV and eo = 3 eV. The band structure 
and Fermi surface at half-filling are shown in Fig. [51^ a) 
and (b), respectively. Below we will take 5sdw = 1-8 eV, 
for which the mean-field equations yield a SDW with 
nesting vector Q = (7r/a,7r/a), critical temperature 
TsDw = 138 K, and T = gap A = 21.3 meV. The sys- 
tem is insulating at T = 0, with the SDW gap completely 
removing the Fermi surface. 

We plot the imaginary parts of the interband, in- 
traband, and total transverse susceptibilties for q = 
{qx,qy = Qx) in Fig. [Tl^a), (b), and (c), respectively. We 
consider first the interband contribution. For q suffi- 
ciently close to Q, we find a continuum of single-particle 
excitations. In contrast to the results for the Hubbard 
model (Fig. ^ , the magnitude of the transverse suscep- 
tibility in this region tends to increase with increasing 
Lu. This can again be explained in terms of the DOS of 
the non-interacting model, which now increases as the 
energy is raised (lowered) away from the Fermi energy 
up (down) to a van Hove singularity at 3eV (— 3eV) 
in the electron-like (hole-like) band. The "density of 
excitations" contributing to the susceptibility therefore 
also increases with u. As the SDW state is insulat- 
ing, with a minimum energy of 2A required to excite 
a quasiparticle across the gap, the continuum is sharply 
bounded at uj — 2 A ~ 42.6 meV. The continuum is also 
bounded by a dispersing V-shaped feature with minimum 
at q = 0.54 Q, which is not seen for the Hubbard model. 
The absence of any weight at small q is anticipated from 
the band structure in Fig. [6l which shows that the mini- 
mum wave vector for an interband transition with energy 
UJ < 400 meV is q w 0.5 Q. The V-shaped feature is plot- 
ted in detail in Fig.[5Ja). As shown in Fig.[5Jb), it is due 
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FIG. 5: Diagrammatic representation of the Dyson equation (|29|l for (a) X-^+'^oo ^'I'i (b) X-'^+'Joo- The curved hues are the 
mean-field Green's functions in the SDW state. When the label b — c, f follows the label a = c, f in the direction of the arrow, 
the corresponding Green's function is G"''. 

with the spin-wave dispersion. 

We now turn our attention to the intraband contribu- 
tion to the transverse susceptibility in Fig. [T^b). Apart 
from a forbidden region close to q = 0, this appears al- 
most like a mirror image of the interband susceptibility, 
albeit much reduced in weight. In particular, we note 
a V-shaped dispersing feature at q = 0.46 Q, the ten- 
dency of Imx'"+'*(q, w) to increase with increasing cu, 
and a dispersing feature at the edge of the q « for- 
bidden region, which resembles the paramagnon close to 
q = Q. The presence of the V-shaped feature is partic- 
ularly interesting, as the discussion above indicates that 
it is due to interband excitations. Thus we find that 
interband excitations give a significant contribution to 
the intraband susceptibility. This is confirmed by exam- 
ining the Dyson equation for x'^+ooi '^f- Tig. [Hb): for 
gcc = 9ff = 52b = <72a = 0, as assumed here, the intra- 
band susceptibilities do not appear on the right-hand side 
of the equation so that the RPA-enhancement of x'!f+ oo 
stems only from the umklapp susceptibilities x^J+^qo ^^'^ 
X'^-+^QQ- The coupling to these terms in the Dyson equa- 
tions is through the anomalous Green's functions G"^^ and 
G^"^, which reflect the mixing of the states in the electron- 
like and hole-like bands separated by the nesting vector 
Q in the SDW phase. Consequently, the intraband sus- 
ceptibility is similar to the interband susceptibility, but 
shifted by Q. 
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FIG. 6: (Color online) (a) Band structure and (b) Fermi sur- 
face of the non-interacting insulating excitonic model. In (b), 
the nesting vector Q = (7r/a,7r/a) is also shown. 



to the weak nesting of the hole band at k = 0.23 Q with 
the electron band at k = 0.77 Q. For the energies con- 
sidered here, to excellent approximation the interband 
susceptibility depends only upon |(5q| = |Q — q|. 

For q near Q, Fig. [7]Ja) shows a spin-wave dispersion 
which appears to intersect the continuum and continue as 
a paramagnon. Figure [9] reveals, however, that the situ- 
ation is more complicated: the spin- wave dispersion does 
not intersect the continuum, but instead flattens out as it 
approaches u = 2A and disappears at q « 0.985 Q. As in 
the Hubbard model, the paramagnon and spin- wave dis- 
persions appear to avoid one another. The paramagnon 
nevertheless seems to connect to significant weight lying 
just inside the continuum region at the intersection point 
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(a) q^a/Tt = qya/TZ 
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(c) = qyCilT^ 



FIG. 7: (Color online) Imaginary part of (a) the interband, (b) the intraband, and (c) the total transverse susceptibility in the 
insulating excitonic model for q = [q^, Qy = qx)- The spin- wave dispersion is visible as the dark line at a; < 2 A close to Q in 
(a) and (c). Note the logarithmic color scale. 



The total transverse susceptibility in Fig. \7lc) clearly 
shows the partial symmetry of the response about q — 
Q/2 but also the asymmetric distribution of weight. 
Imx (j-") for |q| < IQI/2 is roughly one order of mag- 
nitude smaller than at q' = Q — q. 



1. Spin-wave velocity 

The calculation of the spin-wave velocity proceeds as 
for the Hubbard model. For gcc — gjf — 52b = 0, solv- 
ing the Dyson equations for the interband susceptibilties 
in Eq. (|32cp and Eq. (|32d[) again involves the inversion 
of a 2 X 2 matrix, which has the determinant 



V{q,Lo) 



1 



c//c(0) c/c/(0) 
9cfX- + ^00 -ff2aX- + ,00 



1 n ^f'^^'i (") 



52aX- + ,00 



' ^^fcfc{0) fccfiO) 
5c/X- + ,00 +52aX- + ,00 



5c/X- + ,00 +52aX- + ,00 



(37) 



Expanding this determinant about lo — and (5q = Q — 
q = 0, we obtain the low-energy linear form lo — csw \ Sci\ 
of the spin-wave dispersion. For the band structure con- 
sidered here, the spin-wave velocity is given by 



Csw 
where 

ao = 



92a 



{aj + 20002) (g| 



2a252a 



-E 



A2 



1 



Ol 



El 



k 



(38) 



(39) 



(40) 
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FIG. 8: (Color online) (a) Imaginary part of the interband 
transverse susceptibility in the insulating excitonic model for 
q — {qx,qy = qx) close to 0.54 Q. Note the linear color scale, 
(b) Band structure along the Brillouin zone diagonal, showing 
the nesting responsible for the dispersing feature in (a). 



as 



and 



A2 



sin^ k.,a ) (41) 



(42) 



We plot csw s-s a function of .gsDW for different values of 
Qcf in Fig.fTnTa). The behaviour of the spin- wave velocity 
for gsDW ^ i is included as an inset. For psDW > dcf we 
always find Cg-y^ > 0; for sufficiently large g^f > gsBW, 
however, we have Cg^y < which indicates that the sys- 
tem becomes unstable towards a different ground state. 
This is not surprising, as for gcf > 5sdw > > g2a the 
effective coupling 5sdw is smaller than that for the CDW 
or SCDW. In the opposite case g2a > ffSDW > > gcf 
the coupling constants for the SDW and CCDW states 
are equal and always greater than those for the CDW and 
SCDW, and so the SDW remains stable. In Fig. \Mh) 
we plot the spin-wave dispersion in the excitonic model 
as a function of 5q. Compared to a Hubbard model with 
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FIG. 9: (Color online) Imaginary part of the interband 
transverse susceptibility in the insulating excitonic model for 
q ~ (qxyQy ~ Qx) close to Q. The spin- wave dispersion is 
visible as the thick black line in the lower left-hand corner. 
Note the linear color scale. 
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FIG. 10: (Color online) (a) Spin-wave velocity csw as a func- 
tion of (?sDW = gcf + 32a in the insulating excitonic model. 
Inset: csw for a larger range of (?sdw- (b) Spin- wave dis- 
persion in the excitonic model and its low-energy linear form 
for gsDW = 1.8 eV as a function of 5q = q — Q. Shown for 
comparison is the spin-wave dispersion in the Hubbard model 
from Fig. |4jb) for the same T — gap. 
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identical T = gap, the spin-wave dispersion has both 
a higher low-energy velocity and remains approximately 
linear up to higher energies [see Fig. [D^b)] . 

Although Eq. (|38p is a rather complicated function of 
gcf and g2a, for gsDW < 3eV the spin- wave velocity in 
the excitonic model shows remarkably little dependence 
upon the interaction constants, in contrast to the Hub- 
bard model results in Fig.ma). Instead, the value of csw 
is fixed by the band structure: for an excitonic gap A <C t 
(the weak-coupling limit) we have to excellent approxi- 
mation Csw ~ Vi?/\/2 where vp is the average Fermi ve- 
locity. This is anticipated by the results of Refs.[2^and[23l 
for chromium, where it was found that csw = \/veVhJ^ 
where fe(?i) is the electron (hole) Fermi velocity, and the 
factor of l/-\/3 arises because a three-dimensional Fermi 
surface is considered. It is also consistent with our obser- 
vation that X i-(cti'^) is rather insensitive to the choice 

of gcf and 52a for small ^sdw- 

The behaviour of cgw in the strong-coupling regime of 
the Hubbard and excitonic models is also qualitatively 
different. In the former, the U ^ t limit of Eq. ([7]) gives 
Csw = V^J where J = At^ /U is the exchange integral 
of the corresponding effective Heisenberg modelJ^i^i^ 
In the excitonic model, however, the inset of Fig. fTUTal 
reveals that csw ha-s only weak dependence upon the in- 
teraction strength for 5sdw ^ i. A strong-coupling ex- 
pansion of Eq. (j38p gives the limiting value csw = ^/2ta. 
The interpretation of the strong-coupling limit in the ex- 
citonic model is not straightforward: as (?sdw — 00, si- 
multaneous occupation of the c and / states on the same 
site is forbidden, but double occupation of the c and / 
states is allowed. Since we work at half-filling, one might 
expect a checkerboard orbital ordering with filled c states 
on one sublattice and filled / states on the other, which 
is incompatible with SDW order. However, it has been 
shown in a spinless two-band model that such a state is 
unstable towards an excitonic insulator or a phase with 
either the c or / states fully occupied for cq ^ 0^2. How 
this result would change in the presence of spin is not 
clear. In any case, the psDW 3> t limit seems some- 
what unphysical without also considering gcc and 5// to 
be large, and so we do not further discuss the strong- 
coupling regime here. 

The evaluation of the intraband susceptibilities pro- 
ceeds similarly but here the denominator is I?(q-|- Q, oS). 
This yields an identical spin-wave dispersion but shifted 
to q = 0. As in the Hubbard model, however, the spin 
wave has vanishing weight close to the zone center and 
is barely visible as it exits the continuum in the lower 
right-hand corner of Fig. [Tl^b). 



B. Metallic SDW state 

It is more generally the case that the nesting condition 
~ ~^k+Q i^ only approximately satisfied. Further- 
more, there may be portions of the Fermi surface that 
do not participate in the excitonic instability, as is the 
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FIG. 11: (Color online) (a) Band structure and (b) Fermi 
surface of the non-interacting metallic excitonic model. In 
(b), the nesting vectors Q = (7r/a,0) and Q' = (0,7r/a) are 
also shown. 



case in chromium,— i^^i^ The pnictides also have a com- 
plicated Fermi surface involving several bands. Although 
the numerous models for the band structure differ in their 
details r^ii^ii^i^^i^i^i^i^ there is general agreement that 
in the "unfolded" Brillouin zone corresponding to the 2D 
iron sublattice the nesting of hole pockets at k = (0, 0) 
with electron pockets at (tt/o, 0) or (0,7r/a) is primar- 
ily responsible for the SDW. In the physical, tetragonal 
Brillouin zone, both (tt/o, 0) and (0, tt/o) are folded back 
onto the M point, leading to two electron pockets around 
that point. The wave vectors in the present paper re- 
fer to the unfolded zone. Apparently only one of the elec- 
tron pockets undergoes the excitonic instability, yielding 
a SDW with ordering vector Q = (7r/a,Q), say. The 
other electron pocket at Q' = (0, tt/o) remains ungapped. 
We can capture the basic features of this scenario within 
a two-band model by including one hole pocket around 
(0, 0) and treating the two electron pockets as belonging 
to the same band. We thus assume the band structure 



— 2t cos k^a cos ky a ^ 



J - 



2t (cos k^a + cos kyo) + e/. 



(43a) 
(43b) 



We take t = 1 eV, ec = 1.5 eV, e/ — — 3.5eV, and fix 
the doping at n — 1.916, which gives electron and hole 
pockets of identical area. The band structure and Fermi 
surface are illustrated in Fig. [TlTa') and (b) , respectively. 
Note that the hole pocket is nearly but not quite perfectly 
nested with both electron pockets. We impose a single- 
Q SDW with ordering vector Q — (7r/a,0). For gsDW = 
1.873eV we find a mean-field state with critial temper- 
ature TsDW = 132 K and T = gap A = 21.3 meV. In 
the r = SDW state both the hole pocket at the zone 
center and the electron pocket at (tt/a, 0) are completely 
gapped, while the electron pocket at (0, n/a) remains in- 
tact. 

The imaginary parts of the interband, intraband, and 
total transverse susceptibilities for q — [q^^ 0) are shown 
in Fig. [TWa). (b), and (c), respectively. Our results 
are very similar to those for the insulating SDW model 
in Fig. [T] The slightly higher magnitude of the trans- 
verse susceptibility is due to the greater density of states 
in the electron-like band. The similarity is not surprising, 
as the relevant excitations in both models have identical 
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FIG. 12: (Color online) Imaginary part of the (a) interband, (b) intraband, and (c) total transverse susceptibility in the metallic 
excitonic model for q — {q^, 0). The spin- wave dispersion is visible as the dark line at lo < 2 A close to Q in (a) and (c). Note 
the logarithmic color scale. 
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FIG. 13: (Color online) Imaginary part of the interband transverse susceptibility in the metallic excitonic model for q — 
(yr/a — qcosd,qsm.9) for — tv/S. The spin-wave dispersion is visible as the dark line at < 2A close to Q. Note the 
logarithmic color scale. 
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FIG. 14: (Color online) (a) Electron and hole Fermi pockets 
superimposed upon one another, (b) Density of states close to 
the Fermi energy, (c) Comparison of the band structure near 
the intersection of tk+q s^nd in the normal and SDW phases 
for k = (fca;,0). (d) Same as in (c), but for k = {kx,ky = kx). 



origin, i.e., excitations between states close to two Fermi 
pockets which are gapped by an excitonic SDW insta- 
bihty. The states close to the ungapped Fermi pocket 
do not contribute to the intcrband susceptibihty for the 
plotted range of (q, w). Although these states do con- 
tribute to the intraband susceptibility for small values 
of q, they are only negligibly mixed with states in the 
hole-like band, and thus are not RPA-cnhanced by the 
interband interactions. 

In contrast to the insulating SDW state studied 
in Sec. 1111 Al here the intcrband susceptibility docs 
not just depend upon \Sq\: although it is identical for 
q — iqx,0) and (7r/a,7r/a — qx) by tetragonal symme- 
try, and quantitatively very similar along q = (jr/a — 
qx/V^,iqx/V^), away from these high-symmetry lines 
in q-space we find that the continuum can extend to sig- 
nificantly lower energies. This is shown in Fig. 1131 where 
we plot Im X-+''{li ^) for q = {iT/a — q cos 6, qs'm 9) with 
9 = Tr/8. Although the response for oj > 100 meV is very 
similar to that in Fig. [TWa). we see that the lower edge 
of the continuum is not constant at a; = 2 A, but instead 
shows higher and lower thresholds which coincide only at 
special values of q. 

The origin of this threshold behavior is the imperfect 
nesting of the Fermi pockets. Consider Fig. [Tif al. which 



shows the superimposed hole and back-folded electron 
Fermi pockets: along k = {kx,0) or (0, A:^), the width of 
the hole Fermi pocket is greater than that of the elec- 
tron one, while the reverse is true for \i = (kx , ky = kx) 
or (kx,ky = —kx). In the former case, the intersection 
of the non-interacting electron-like and hole-like disper- 
sions therefore occurs above the Fermi energy [Fig. [T4r c)]. 
while in the latter it occurs below the Fermi energy 
[Fig. [MTd)]. In the reconstructed bands of the exci- 
tonic model, the SDW gap is always centered at the 
point of intersection of the original bands, as can be seen 
from Eq. ([^5]) and in Figs. fHf c) and (d). In general, 
the difference between the Fermi energy and the bot- 
tom of the reconstructed electron-like band, A+ k, and 
the difference between the Fermi energy and the top of 
the reconstructed hole-like band, A_.k, will be unequal. 
The minimum energy for an interband excitation with 
wave vector Q + ^q is therefore mink(A±.k -I- Azp.k+iq)- 
For Jq along the high-symmetry directions mentioned 
above, the tetragonal symmetry of the Fermi pockets 
ensures that this minimum energy is 2 A. Away from 
these directions, however, the energy difference depends 
upon (5q. For example, the states near the Fermi sur- 
face in Fig. [MTa'l at k = {0,ky) and {kx,ky = kx) 
are connected by Sq = (g cos(7r/8), — gsin(7r/8)) with 
q = 0.17 7r/a; from Fig. [MF c) and (d) we see that the 
minimum energy for single-particle excitations with this 
wave vector is 1.2 A = 25.6 meV, which marks the lowest 
edge of the continuum in Fig. 1131 The upper threshold 
originates from the maximum energy connecting the top 
of the hole-like band and the bottom of the electron-like 
band, which for this q is 2.8 A = 59.6 meV. For the re- 
mainder of this paper we shall restrict ourselves to high- 
symmetry directions. 

Since only one electron pocket is gapped, the direc- 
tions (qxiO) and {0,qy) are not equivalent. The imagi- 
nary part of the interband transverse susceptibility along 
q = (0, qy) is shown in Fig. [151 This is the direction to- 
wards the second nesting vector Q' = (0, tt/o), which was 
not selected by the SDW instability. For q sufficiently 
close to Q', we thus find the response generated by transi- 
tions between states near the (gapped) hole Fermi pocket 
and states near the ungapped electron Fermi pocket. At 
Lu > 2 A, this is very similar to the interband suscepti- 
bility near Q [Fig. [T2r a)] , reflecting the small changes to 
the band structure at high energies upon opening of the 
SDW gap. The differences are more striking at lower en- 
ergies. In particular, comparing Fig.[T5]to Fig.fT^a). we 
see that the continuum extends to lower energies close to 
Q' than close to Q. The minimum energy required for 
a single-particle excitation between the states near the 
ungapped electron pocket and the gapped hole pocket is 
smaller than 2A, thus giving a lower threshold for the 
continuum near Q'. 

Another significant difference concerns the spin-wave 
dispersion. The spin- wave dispersion near Q is visible in 
the lower left-hand corner of Fig. [T^ a). and it intersects 
the continuum and appears to continue as a paramagnon. 
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FIG. 15: (Color online) Imaginary part of the interband transverse susceptibility in the metallic excitonic model for q = (0, qy) 
The spin- wave dispersion is visible as the dark line at a; < 2 A close to Q' {qy = 0). Note the logarithmic color scale. 
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FIG. 17: (Color online) Spin-wave dispersion in the metallic 
excitonic model close to Q for 5q — {—qcos9,qsm0) with 
9 — (solid red line) and 9 — n/8 (thin dotted red line). 
The dispersion for 9 = 7r/4 is indistinguishable from the 9 = 
case. We also plot the spin-wave dispersion close to Q' 
(dashed black line). Inset: imaginary part of D(q, a;) for each 
dispersion. We have used a 20000x20000 k-point mesh and a 
width S = 0.1 meV to calculate the mean-field susceptibilities. 
Note that the finite value of Im©(q, w) for the spin waves 
close to Q for uj < 2A is an artifact of the finite width S. 



FIG. 16: (Color online) Imaginary part of the interband trans- 
verse susceptibility in the metallic excitonic model for (a) 
q = {qx,0) close to Q and (b) q = {0,qy) close to Q'. In 
both panels the spin-wave dispersion is visible as the thick 
black line in the bottom left-hand corner. Note that in (b) 
that the continuum region starts at a; ~ 0.6 A. In both panels 
we use a linear color scale. 



From Fig. [13 we see that there is also a gapless Gold- 
stone mode at q = Q'. This mode is gapless since it 
rotates the single-Q SDW into a superposition of Q and 
Q' SDWs, which is degenerate with the single-Q SDW 
in our tetragonal model. Although there appears to be 
a spin- wave branch around Q', it is not as distinct as 
in Fig. [TWa) due to the lower threshold of the continuum. 



We therefore plot the interband transverse susceptbilities 
for a finer q-resolution near Q and Q' in Fig. [TBFa) and 
(b), respectively. As expected from the discussion above, 
the former is qualitatively identical to Fig. [51 The latter, 
in contrast, shows several novel features: the spin- wave 
dispersion does not curve away from the edge of the con- 
tinuum but rather intersects it with little change in veloc- 
ity and the spin-wave and paramagnon features approach 
much closer to one another than for q w Q. Although 
it is not clear from Fig. riHT b). the spin- wave and para- 
magnon dispersions do not intersect, and the spin waves 
become damped at uj fa 1.7A. 

To obtain the spin-wave dispersion, we must again 
solve Rel?(q,w) = with I?(q,cj) given by Eq. (|57)l . 
We have not been able to obtain analytical expressions 
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FIG. 18: (Color online) Imaginary part of the transverse sus- 
ceptibility as a function of uj at various values of q near (a) 
Q and (b) Q' in the metallic excitonic model. We have cal- 
culated the mean-field susceptibilities using a 30000 x 30000 
k-point mesh and a width 5 = 0.2 meV. 



qy = 0.99 7r/a. To examine this more closely, we plot 

Imx h(Qi^) a function of uj for fixed q near Q and 

Q' in Fig.fTWa) and (b), respectively. In the former case, 
we see the step-like start of the continuum aX oj — 2 A. 
The peak at this energy is due both to the remnant of 
the spin-wave branch (at least for qx — 0.98 7r/a) and to 
the enhancement of the DOS at the edge of the SDW 

gap. The finite value of Imx h(lj'^) fo'" ^ < 2A is an 

artifact of the finite width 5. The susceptibility near Q' 
is qualitatively different: the lower threshold of the con- 
tinuum is at u>i = 0.65 A, and immediately above this 
the susceptibility increases continuously as ^oj — oji. At 
llI2 — 1.3 A, the susceptibilility abruptly starts to increase 
more steeply. The locations of these two thresholds cor- 
respond to the kinks in the spin-wave dispersion. The 

rapid increase of Imx hCli ^) above uj2 accounts for the 

strong increase in the damping (see inset of Fig. I17|) . 

As for the interband susceptibility near Q, the ori- 
gin of the lower threshold is the variation of A±^k- The 
difference is that here the threshold originates from the 
minimum energy required for a single-particle excita- 
tion between the states near the gapped hole pocket 
and the states near the ungapped electron pocket, uji = 
mink(A+_k, A__k)- From Fig. fWc) and (d) we de- 
duce LOi ~ 0.6 A, closely matching the lower threshold 

in Fig. fTSF b). The strong increase in Imx [-(Qi ^) above 

L02 is due to the peaks in the DOS located at ±ci;2 on 
either side of the Fermi energy, shown in Fig. [MFb) : be- 
cause of this DOS enhancement, the "density of excita- 
tions" between states close to the gapped and ungapped 
Fermi pockets is increased above W2- 



for the spin-wave velocity, however, as the Fermi distribu- 
tion functions appearing in the mean-field susceptibilities 
cannot be expanded as a Taylor series in 5(\ due to the un- 
gapped electron Fermi pocket. Plotting the dispersions 
at q « Q and q « Q' in Fig. [ITl we see that the velocity 
at Q is roughly 25% higher than at Q'. Despite the vari- 
ation in A-t^k, there is no anisotropy of the low-energy 
spin- wave velocity. The difference between the results for 
9 = {] and = 7r/8 at higher energies is due to the lower 
edge of the continuum in the latter case. Whereas the 
spin waves close to Q have a very similar dispersion com- 
pared to the insulating model [Fig.fTCTb)]. the dispersion 
close to Q' has two noticeable kinks at oj = 0.65 A and 
w = 1.25 A. As shown in the inset, these kinks coincide 
with abrupt changes in Iml?(q, cj): Iml?(q, w) becomes 
finite at a; = 0.65 A, and starts to sharply increase at 
OJ = 1.25 A. The first feature corresponds to the onset 
of Landau damping as the spin-wave branch enters the 
continuum. The second feature is a result of the DOS, 
as discussed in the following paragraph. 

Examining the lower edge of the continuum in both 
panels of Fig. [TBI we see that whereas the continuum dis- 
appears sharply at w = 2 A near Q, it appears to vanish 
more smoothly near Q'. In the latter case there are two 
distinct thresholds, which are particularly visible around 



IV. EXPERIMENTAL SITUATION 

This work ultimately aims to shed light upon the na- 
ture of the antiferromagnetism in the iron pnictides, in 
particular the extent to which it is itinerant or local- 
ized in character. There are several published results of 
inelastic neutron scattering examining the spin excita- 
tions in the antiferromagnetic state of CaFe2 As2 ^"^-I i^'^i^^ 
SrFe2As2r^ and BaFe2As2 4^'^ In these experiments, 
only transverse excitations contribute to the neutron- 
scattering cross section, allowing us to write it as 

(X |i^(q)p [nB{uj) + 1] Imx-+(q,c.), (44) 

where F(q) is a form factor and ris(aj) is the Bose- 
Einstein distribution function. A direct, quantitative 
comparison between theory and experiment would re- 
quire a more realistic model for the low-energy band 
structure than the one we are using. We neverthe- 
less make several general remarks relating what we have 
learnt about the spin excitations in the excitonic SDW 
state to the experimental results. 

We first review the experimental situation. Despite 
considerable variation in the Neel temperature within 
the y4Fe2As2 {A =Ca, Sr, Ba) family, the static magnetic 
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properties of these compounds are rather similar. In par- 
ticular, antiferromagnetism only occurs in the presence 
of an orthorhombic distortion, which fixes the ordering 
vector Q. Experiments on the low-energy spin dynamics 
are also in broad agreement: there is a strongly dispering 
spin wave close to Q ^46,47,54,55,56,57 ^jjg gpin-wave velocity 
is anisotropic jili^i^i^i^ and the spin-wave dispersion 
has a gap of energy 6-10 meV4^'^ i^'^i^^ '^i^ At present, 
however, there is considerable disagreement over the 
high-energy excitations. For CaFe2As2, it was reported^ 
that the spin wave is strongly damped at energies above 
100 meV, suggesting the presence of a particle- hole con- 
tinuum. On the other hand, although Zhao et alM found 
similar spin- wave velocities, they did not observe any sig- 
nificant jump in the damping of the spin wave below 
200 meV, which would indicate the intersection of the 
spin wave dispersion with the continuum. The results 
for BaFe2As2 show greater inconsistency, with reports'*^ 
of strong spin excitations possibly up to 170 meV in stark 
disagreement with claims of spin-wave damping by con- 
tinuum excitations at energies as low as 24 meV.— 

The results of Refs. [sj and are most consistent 
with itinerant antiferromagnetism, as the existence of a 
particle-hole continuum is a key feature of this scenario. 
Interpreting the latter experiment^ in terms of the cxci- 
tonic model, we deduce a SDW gap of A « 50 meV. This 
is nearly twice the estimate A « 30meV of the T — 
gap based on ARPES for SrFe2As2 3 Although a SDW 
gap of only 12meV for BaFe2As2, which we could infer 
from Ref. [53, seems low, we have seen above that spin- 
wave damping sets in at energies much smaller than 2 A, 
depending upon the details of the reconstructed band 
structure. In order to fit the results of Ref. i47i into the 
excitonic picture, however, we require a SDW gap of at 
least 100 meV implying a rather high value of the ra- 
tio A/fcsTsDw ^ 7. These results instead support a 
local-moment picturc i^^i'^^'^^i'*^ The absence of the con- 
tinuum is nevertheless surprising since ARPES shows 
clear evidence for quasiparticle bands at low energies, 
which suggests a possible resolution i^°i^^ the imperfect 
nesting of the elliptical electron pockets with the circu- 
lar hole pocket is expected to yield incompletely-gapped 
Fermi surfaces in the SDW state, which implies that con- 
tinuum excitations are present down to zero energy. As 
such, the spin waves would be damped at all energies, 
and the jump in the damping characteristic of the en- 
try into the continuum is absent. Such an explanation 
is of course at odds with Ref. IH, indicating the need for 
further work to clarify the experimental situation. 

The reported 40% anisotropy of the spin- wave velocity 
within the ah plane^i^ is quite remarkable. Although 
this effect is absent from our results due to the tetrago- 
nal symmetry of the Fermi pockets, it nevertheless seems 
rather too large to be accounted for by the expected el- 
liptical shape of the electron-like Fermi pockets in the 
pnictides42. Experimental results also do not show a sec- 
ond spin- wave branch at Q', as found here for the metal- 
lic SDW model. Both observations are likely due to the 



orthorhombic distortion in the SDW phase, which lifts 
the degeneracy of the (tt, 0) and (0, tt) SDWf2£ and do 
not imply a failure of the excitonic scenario. 

We finally remark upon the gap in the spin-wave dis- 
persion in the pnictides. Due to the absence of of mag- 
netic anisotropy is our model, we always find Goldstone 
modes in the SDW phase. As demonstrated in Fishman 
and Liu's study of manganese alloys,.^^ a gap is possible in 
an excitonic SDW state in the presence of magnetoelastic 
coupling. The magnetoelastic coupling in the pnictides is 
indeed strong, as evidenced by the role of the orthorhom- 
bic distortion in fixing the polarization and the ordering 
vector of the SDW^^iii^ suggesting that it might be re- 
sponsible for the spin- wave gap. 

In summary, the neutron-scattering data for the anti- 
ferromagnetic state in the pnictides are currently unable 
to decide upon the origin and character of the magnetism. 
We have shown that the excitonic SDW scenario gives 
spin-wave excitations in qualitative agreement with ex- 
periments. An obvious direction of future work is there- 
fore to examine the spin excitations based on more real- 
istic band structures. Considering the imperfect nesting 
of the electron and hole pockets in the pnictides, it will 
be particularly interesting to address the possibility of 
incommensurate SDW order.— The effects of the inter- 
actions not directly contributing to the SDW instability 
should be included. Comparison of our results with those 
obtained within a model explicitly accounting for the or- 
bital character of the bands is also important. Further- 
more, the orthorhombic distortion and a magnetoelastic 
coupling should be implemented for greater realism. Al- 
though the spin excitations in more sophisticated models 
will differ in their details from those presented here, we 
nevertheless think that our results will remain qualita- 
tively correct and will thus be valuable in interpreting 
future experiments. 



V. SUMMARY 

We have presented an analysis of the zero-temperature 
transverse spin excitations in the excitonic SDW state 
of two-band, 2D models with nested electron-like and 
hole-like Fermi pockets. Using the RPA, we have derived 
the Dyson equation for the spin susceptibility and have 
shown that the total spin susceptibility can be divided 
into contributions from interband and intraband excita- 
tions. We have solved the Dyson equation in the special 
case when only the interactions responsible for the SDW 
are non-zero. While the interband excitations are then 
directly enhanced by the interactions, the intraband exci- 
tations are still indirectly enhanced due to the mixing of 
the electron-like and hole-like states in the SDW phase. 
The susceptibility exhibits collective spin- wave branches 
close to the SDW ordering vector Q and also, with much 
smaller weight, close to q = 0, as well as a continuum of 
single-particle excitations at energies above a threshold 
of the order of the SDW gap. 
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Depending upon the non-interacting band structure, 
the opening of the excitonic gap can result in quahta- 
tively different SDW states. This has been illustrated 
by considering two models, one which becomes insulat- 
ing in the SDW state and another which remains metal- 
lic due to the presence of an ungapped portion of the 
Fermi surface. For comparison, we have also performed 
the corresponding calculations for a 2D Hubbard model 
with the same mean-field SDW gap. Differences in the 
spin excitations between the insulating and metallic mod- 



els occur only at low energies and mainly close to the 
nesting vector Q' between the (gapped) hole pocket and 
the ungapped electron pocket, which is essentially unaf- 
fected by the SDW formation. We have also discussed 
data from neutron-scattering experiments in light of our 
results. We conclude that the available data do not yet 
allow us to distinguish between an excitonic SDW and a 
local-moment scenario for the antiferromagnetic order in 
the pnictides. 
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